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ABSTRACT 

We report the discovery of a bright (V ~ 13.7), extremely low-mass white dwarf in a close double 
degenerate system. We originally selected GALEX J171708. 5+675712 for spectroscopic follow-up 
among a group of white dwarf candidates in an ultraviolet-optical reduced proper-motion diagram. 
The new white dwarf has a mass of 0. 18 Mq and is the primary component of a close double degenerate 
system (P = 0.246137 d, K\ — 288 kms" 1 ) comprising a fainter white dwarf secondary with Mi ps 
0.9 Mq. Light curves phased with the orbital cphemeris show evidence of relativistic beaming and 
weaker ellipsoidal variations. The light curves also reveal secondary eclipses (depth w 8 mmag) while 
the primary eclipses appear partially compensated by the secondary gravitational deflection and are 
below detection limits. Photospheric abundance measurements show a nearly solar composition of 
Si, Ca, and Fe (0.1 — 1 0), while the normal kinematics suggest a relatively recent formation history. 
Close binary evolutionary scenarios suggest that extremely low mass white dwarfs form via a common- 
envelope phase and possible Roche-lobe overflow. 

Subject headings: binaries: close — stars: individual (GALEX J171708. 5+675712) — white dwarfs 



1. INTRODUCTION 

Many new members of the class of extremely low- 
mass (ELM) white dwarfs were discovered in the 
Sloan Digital Sky Survey (SDSS) and the Luyten 
surveys (e.g., Luyten Palomar, Luyten Half-Second, 
New Luyten Two-T enths), such a s SDSS J123 4-0228 



(jLiebert et alj l200l, LP400-200 (iKawka et alj [2006), 
SDSS J0917+4638 kilic eTail 120071) . and NLTT 11748 



(|Kawka k. Venned I2009D . These objects are the prod- 
ucts of close binary evolution and likely emerged from 
a co mmon-envelope phase a nd possible Roche-lobe over- 
flow (jNelemans et al.l 120011 ). Some ELM white dwarfs 
are paired with a neutron star companion, e.g., the pul- 
sar J1012+5307 (|van Kerkwiik et all [1990) . implicating 
a range of initial secondary masses extending to > 8 Mq . 

We have identified ~ 200 bright hot subdwarf and 
white dwarf stars in a joint Galaxy Evolution Ex- 
plorer {GALEX) ultraviolet (UV) surve y and Guide 
Star Catalogue (GSC) survey (paper I: I Vennes et alJ 
l20lTl and paper II: Nemeth et al., in preparation). 
A reduced proper-motion diagram based on the pho- 
tographic color V ~ Vgsc (Hv = V + 51og/i + 5) 
versus the color index A+jv — Vgsc segregated the 
hot subdwarf and white dwarf candidates. As part 
of a second installment of our catalog of UV-selected 
stars, we identified GALEX J171708.5+675712 (here- 
after, GALEX J1717+6757) as a likely white dwarf with 
H v « 12.3. 

We present photometric and spectroscopic observa- 
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tions (Section 2), and a model atmosphere and peri- 
odogram analysis showing GALEX J1717+6757 to be a 
new ELM white dwarf in a close binary and we analyze 
its photometric properties (Section 3). In Section 4, we 
estimate the distance, kinematics, and age of the system 
and discuss some implications for close binary evolution- 
ary scenarios. 

2. OBSERVATIONS 

We observed GALEX J1717+6757 on UT 2011 Jan- 
uary 28 using the Ritchey-Chretien Focus Spectrograph 
(RC-spec) attached to the 4-m telescope at Kitt Peak 
National Observatory (KPNO). We employed the T2KA 
CCD and the KPC-10A grating (316 lines per mm) with 
a dispersion of 2.75 A per pixel in first order and cen- 
tered on 5300 A. We also inserted the order-sorting filter 
WG360 and covered a useful spectral range from 3600 
to 7200 A. The slit width was set at 1.5" resulting in 
a spectral resolution of ss 5.5 A. We obtained two con- 
secutive exposures of 420 s each. A preliminary model 
atmosphere analysis revealed a relatively hot and sublu- 
minous object in a region of the HR diagram populated 
with ELM white dwarfs. We also measured a radial ve- 
locity ~ 300 km s" 1 indicating either peculiar kinematics 
or an orbital motion characterized by a large velocity 
amplitude. 

We therefore obtained time series spectroscopy from 
2011 March 18 to 21, using the 2. 4-m Hiltner telescope 
and modular spectrograph (modspec) at MDM Observa- 
tory on Kitt Peak. A 1024 2 SITe CCD gave 2 A pixel" 1 
and 3.5 A resolution from 4660 to 6730 A. We obtained 
a total of 65 exposures of 480 to 600 s each, spanning 
nearly 6 h of hour angle and covering all phases of the 
orbit, and measured radial velocities of the Ha and H/3 
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TABLE 1 
Photometry and astrometry. 



Parameter 


Measurement 


GALEX F vv 


13.44±0.10 mag 


GALEX Af uv 


13.80±0.10 mag 


SDSS u 


13.548±0.003 mag 


SDSS g 


13.432±0.006 mag 


SDSS r 


13.441±0.003 mag 


SDSS i 


13.701±0.003 mag 


SDSS * 


13.987±0.004 mag 


2MASS J 


13.61±0.02 mag 


2MASS H 


13.69±0.03 mag 


2MASS K 


13.64±0.05 mag 


R.A. (J2000) 


17 17 08.865 


Dec. (J2000) 


+67 57 11.33 


fi a cos 6 


-6.7 ± 2.2 masyr -1 




51.6 ± 0.5 masyr - 1 



lines using a convolution algorithm (Schneide r fc Yound 
1980). We also obtained eight velocities using the same 
setup 2011 June 21 through 25 to further constrain the 
ephemeris (Section 3.2). 

Next, we obtained Johnson BVI photometric time se- 
ries on 2011 February 8 and 21, and March 19 and 22, 
using the 0.6-m telescope at Bialkow Observatory. The 
exposure times were 60 s in V and /, and 100 s in B. 
The data were corrected for second-order extinction ef- 
fects by fitting a line to the differential magnitude versus 
air-mass diagram and subsequently removing the linear 
dependence. Each time series provided approximately 
5 hours of continuous coverage in the three photomet- 
ric bands. We also obtained high signal-to-noise ratio 
white-light time series on 2011 June 2, 3, 4, 6, and 7 
with exposure times of 30 s. The time series provided 
between 2 to 5 hours of continuous coverage. 

Finally, we collected photomet ric measurements from 
the SDSS (lAbazaiian et alj|2009l). 2 Micron All-Sky Sur- 
vey (2MASS, iSkrutskie et all 120061) and GALEX sur- 
vey (Table [J). We corrected the GAL EX photometry 
for non-linearity (Morrisse v et al.l I2007D and estimated 
an error of 0.1 mag in the process. These measure- 
ments are useful in building a spectral energy distribu- 
tion (SED) and in constraining the atmospheric param- 
eters of the white dwarf (Section 3.1). The SDSS (/-band 
photometric measurement shows discrepancy with other 
photometric measurements and was rejected from the 
SED. Table Q] lists astrometric measurements from the 
Third U.S. Naval Obs ervatory CCD Astrograph Catalog 
(|Zacharias et al.|[2010[ ) that will be used to determine the 
kinematical properties of GALEX J1717+6757 (Section 
4). 

3. ANALYSIS 

3.1. Properties of the primary 

We computed a grid of models in non-local thermody- 
namic equilibrium (non-LTE) using Tlusty- Synspec 
(|Hubenv fc Land Il995t ILanz fc Hubenvl Il995l) versions 
200 and 48, respectively. When fitting the Balmer lines 
(Ha to H13), the best-fitting parameters 

T cff = 14 900 ± 200 K, log g = 5.67 ± 0.05, 

and a low helium abundance log [n(He)/n(H)] < —2.3 
(< 0.05 0, where is the solar abundance relative to 



hydrogen) confirmed the nature of GALEX J1717+6757 
as an ELM white dwarf. Using these parameters we 
estimated the absolute magnitude in the Sloan r band 
M r = 6.99 ± 0.08, and with the corresponding distance 
modulus r-M r = 6.45±0.08 (r = 13.441) we estimated a 
distance of 195 ±7 pc. Figure [T] shows the SED and best- 
fitting non-LTE model; GALEX J1717+6757 lies toward 
a low-density line- of-sight in the inter stellar medium 
(E B -v ~ 0.03 mag. iSchlegel et alj[l998l) with negligible 
extinction in the intervening medium. The mass and age 
of the white dwarf may be estimated by lo cating it on the 
evolu tionary tracks with solar metallichy (Serenclli et al. 
2001 1 ) or with low metallicity (Z = 0.001. jSerenelli et al. 
2002). The mass is model-dependent and the error is 
dominated by the choice of the model grid. Using low- 
metallicity models we find Mi = 0.184 ± 0.001M Q , but 
the mass derived from solar-metallicity models is less 
certain because of a large gap between the 0.169 and 
0.196M Q sequences. Extrapolating from the sequence of 
low mass models (0.148, 0.160, 0.169M Q ) the mass is esti- 
mated at O.185M0. Based on these two model sequences 
and considering that the ELM lies within the available 
sequences, i.e., 0.169 and 0.196M Q , we conservatively es- 
timate Mi = 0.185 ± 0.010M Q and Rx = 0.10 ± O.Oli?© 
for both low- and high-metallicity models. We examine 
the age problem in Section 4. 

The surface composition of GALEX J1717+6757 is 
that of a heavily polluted white dwarf. Figure [T] (in- 
set) shows photospheric heavy element lines that are 
suitable for abundance measurements. Formally clas- 
sified as DAZ white dwarfs, the hydrogen- rich white 
dwarfs with heavy-element spectral li nes constitute ap- 
proximately 20% of the DA population (|Zuckerman et all 
12001 [Koester et all [2005] ). We fitted the Ca IIA3933, 
Fe IIAA4923, 5018, 5169, and Si IIAA5041,5056,6347,6371 
A lines with high-metallicity non-LTE models computed 
using Tlusty-Synspec and measured: 

log [n(Si)/n(H)] = -5.5 ± 0.3 (0.1 0), 
log[n(Ca)/n(H)] = -5.7 ±0.3 (1.0 0), 
log [n(Pe)/n(H)] = -5.1 ±0.3 (0.3 0). 

The atmospheric composition is somewhat below the 
solar composition, with a helium abundance less than 
5% of the solar abundance. High-dispersion and high 
signal-to-noise ratio spectra are required to constrain the 
abundance of other heavy elements. The origin of these 
heavy elements may be fossil or accreted from the im- 
mediate circumstellar en vironment (see a discussion by 
iKoester fc Wilkenll2006h . A direct link to the metallic- 
ity of the progenitor cannot be established because the 
surface composition of white dwarf stars is determined 
by diffusion processes characterized by time-scales much 
shorter than evolutionary time-scales. 

3.2. Binary properties and nature of the secondary 

A periodogram constructed from the radial velocities, 
selected an unambiguous period, and a sinusoidal fit of 
the form v(t) = 7 + K sm[2n(t - T )/P] gave 

T = HJD 2455641.9309 ± 0.0004, 
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Fig. 1.— Spectral energy distribution (2MASS JHK, SDSS ugriz, GALEX N uv and F uv ) of GALEX J1717+6757 compared to the 
best-fitting model 14 900 K and logcj = 5.67. The effect of minimal extinction (Eg_y = 0.005 mag) is added to the model. In the inset, 
the co-added KPNO spectrum shows the upper Balmer lines and a weak photospheric Ca IIA3933 line, and the co-added MDM spectrum 
shows Si IIAA5041, 5056, 5186 and Fe IIAA4923, 5018, 5169 lines. 



P = 0.246137 ±0.000003 d, 

K = 288 ±2 km s~\ 



7 : 



-21.3 ± 1.7 km s~ 



with an RMS residual of only 8 km s -1 . The initial epoch 
To corresponds to the inferior conjunction of the primary. 

Figure [2] shows the spectroscopic trail obtained at 
MDM Observatory: Ha and Si IIAA6347, 6371 clearly 
trace the primary orbit bolstering the photospheric iden- 
tification of the silicon lines. We found no spectroscopic 
signatures of the faint secondary star. The folded radial 
velocity curve follows a circular orbit (Fig. [2]). 

The period and if-velocity imply a mass function for 
the secondary 



/(Ma 



0.609 ± 0.012 M, 



or, assuming Mi = 0.185 ± 0.010 Af Q , M 2 > 0.86 M , 
so the companion must be a compact object. If the 
unseen companion is below the Chandrasekhar limit 
M 2 < 1.35 M Q , then the inclination i > 60°. No ra- 
dio sources are found in the vicinity of the ELM white 
dwarf in the NRAO-VLA Sky Surve y (NVSS) with a sen - 
sitivity of ~ 0.5 mJy at 1400 MHz (|Condon et al.|[l998h . 
contrasting with the detection of t he milli-second pul- 
sar PSR J1012+5307 in the NVSS (Kaplan et alJ [l99cl 
and other radio observations ([Nicastro et al.lll995l ) . Pul- 
sar companions to E LM white dwarfs are relatively rare 
(|Agueros e t al. 2009]) and we conclude that the compan- 



ion to the ELM white dwarf GALEX J1717+6757 is a 
massive white dwarf. 

We folded the Bialkow photometric time series on our 
ephemeris. Including all V data, the time series folded on 
half the period has a weak amplitude of 1.5 ± 0.4 mmag 
and a maximum at phase 0.30 ± 0.02 (statistical error 
only), i.e. close to phase 0.25 expected for ellipticity 
effect. Using "JKTEBOP"0 (IPopper fe Etzell[l98l we 
verified that the amplitude calculated from the model 
for the system parameters agrees with the observed am- 
plitude. The V time series folded on the period has a 
stronger amplitude of 2.2±0.4 mmag with a maximum at 
phase 0.81±0.03, i.e., also « 2a from phase 0.75 expected 
for relativistic beaming. Figure [3] shows the folded V 
time-series and a model curve combining both effects: we 
added stellar distortion with an amplitude of 1.5 mmag 
and relativistic beaming w ith a n amplitude of 2.2 m mag. 
Ivan Kerkwiik et all (|2010D and lShporer et all (|2010D dis- 
cuss the beaming effect in close binaries; because the 
primary outshines the secondary, the amplitude of the 
effect can be written simply as A « a e a /(e a — 1) x (K/c), 
where c is the speed of light, K {= 288 km s -1 ) the ve- 
locity semi- amplitude, and a — hv/kT c g where v is the 
frequency of the photometric bandpass and T e g is the 
stellar effective temperature. In the V band, the am- 
plitude of the effect in GALEX J1717+6757 is expected 
to be A — 2.3 mmag. This is the second star in which 



5 Available at http : //www . astro . keel e ■ ac ■ uk7"| " j kt/codes/ 
jktebop.html 
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Fig. 2. — (Left) trailed spectra showing Si IIAA6347, 6371 and Ha lines tracing the orbital motion (Section 3.2) of the bright ELM white 
dwarf primary. The upper panel is scaled to bring out the Si II features, and the lower to show the core of the Ho line. (Right) radial 
velocity measurements phased on the orbital period with the best-fitting sinusoidal curve tracing the primary orbital motion. The data are 
repeated for one cycle for continuity. The full triangles near phase 0.2 are from the Kitt Peak 4 m, and all other velocities are from the 
MDM 2.4 m. 



relativistic be aming was detected from the ground after 
NLTT 11748 (|Shporer et alllMtOft . 

Adopting for the primary and secondary Mi = 
0.185M©, Ri = O.LRq (Section 3.1) and M 2 = 
0.9Mq, i?2 = O-OO&Rq (us ing mass-radius relations of 
IBenvenuto fc Althauslll999l ). respectively, the secondary 
eclipse would last up to 400 s at i = 90° with a depth of 
7 mmag assuming an identical surface flux for the pri- 
mary and secondary stars. Employing only the June 
data (white light) we observed four secondary eclipses 
at T(HJD 2455000+) = 715.5252, 716.5098, 717.4944, 
and 720.4476, and lasting on average 80-150 s (i = 
86.75 ± 0.05°) with a depth of 8 mmag. The eclipse 
times are consistent with phase = with an offset of only 
—0.003 ± 0.004, and suggest only a small adjustment to 
the orbital period of —0.0000025 d, i.e., within quoted 
errors. Figure [3] shows the white-light time series phased 
with the adjusted period. The primary eclipse, i.e., the 
transit of the secondary, is not clearly detected with a 
maximum depth of w 4 mmag. The primary eclipse is 
partially compensated by gravitational deflection by the 
massive white dwarf. With the adopted parameters and 
a binary separation of 1.7i?0, the Einstein radius of the 
secondary is Re = 0.0036-R©. The total tra nsit depth in- 
cluding the lensing effect may be written as (|Marshll2001[ ) 
A?n w (R 2 /Ri) 2 -2(R E /Ri) 2 : where "1" designates the 
ELM primary, and "2" the massive unseen companion. 
Consequently, the predicted transit is Am w 4 mmag, 
consistent with the deepest feature found at phase 0.5. 



This lensing effect h as been observed in NLTT 11748 
(jSteinfadt et al.|[2010h . 

4. DISCUSSION 

The ELM GALEX J1717+6757 constitutes an in- 
teresting test of evolutionary models. Following low- 
metallicity evolutionary tracks, the white dwarf has a 
cooling age of 88O+3Q Myr, but it is not possible to 
derive an age using the solar-metallicity tracks due to 
a large gap in properties above 0.169M©. Figure [4] 
compares the evolutionary models to samples of ELM 
white dwarfs. We selected the old systems LP 400-2 2 
(jKawka et al.ll2006HKilic et aJJ2009H Vennes et al.ll2 0091 , 
NLT T 11748 dKawka fc Vennesl l2009t iSteinfadt et al.l 
120101 ; [Kawka et al. 2010) and a sample of low-mass white 
dwarfs in the SPSS su rvey (Kulka rni fc van Kerkwiikl 
l2010HBrown et al.|[20Tot iKilic et al 1 1 201 lallbf) . The ELM 
precursor HD 188112 (jHeber et al.l 120030 is also shown 
with a pre- flash model describing low-mass subdwarf 
stars (|Driebe et al.lll99ll ). The evolutionary tracks cover 
the population with masses ranging from > O.15M to 
w O.45M . The ELM white dwarfs (M < 0.2M Q ) oc- 
cupy a relatively high-luminosity region of the d i agram 
characterized by residual burning. lAlthaus et al.l (2001) 
estimates that, taking into account the effect of diffu- 
sion, objects with masses > 0.18M© will experience a 
thermonuclear flash that delays their cooling. The ef- 
fect is also metallicity dependent, with higher-metallicity 
objects experiencing a flash at a comparatively lower 
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Fig. 3. — (Top panel) white-light time series (30 s bins) folded on the cphcmeris (see text) covering the primary (phase = 0.5) and 
secondary (phase = 0) eclipses, and (bottom panel) Johnson V time series (open squares) and model curve (thick line). The upper curves 
in the bottom panel are offset by —30 mmag and depict ellipsoidal effect at half-period and relativistic beaming at full-period (see text), 
while the lower curves are offset by 30 mmag and represent the binned (0.05 phase) light-curve (histogram with error bars) and the model 
curve (thick line). 



mass. Independent age and mass estimates are possi- 
ble in bina ries comprising an E LM white dwarf and a 
pulsar (see lAlthaus et al.l 12001) but arc more difficult 
to obtain in double white dwarf systems. The ELM 
GALEX J1717+6757 is one of the youngest known class 
members, but the estimated age of the system remains 
model-dependent with the metallicity of the progenitors 
still unconstrained. 

The ELM primary in GALEX J1717+6757, with an 
estimated radius R\ — 0.1 Rq fills in only abo ut 1.2% of 
its R oche volume (i?i,L = 0.43 i?©, following lEggfetonl 
1983) resulting in low-amplitude ellipsoidal variations 
compounded with a dominant relativistic beaming effect. 
The system is expected to merge (t melge ) in 7.0 Gyr as - 
suming a secondary mass M 2 = 0.9 M Q (IRitteri ri986K 
and the product of the merger will likely be an ultra- 
massive white dwarf with M w 1 . 1 M© . 

Using a distance d — 195 ± 7 pc (Section 3.1), 
and the velocity and proper motion measurements 
we computed the (U , V, W) velocity vector following 
I Johnson fe Soderbloml (|1987l) : 

(U, V, W) = (-35 ± 3, -21 ± 3,-8 ± 3) kms" 1 . 
The velocity components are characteristic of the thin 



Galactic disc a nd similar to the white dwarf po pula- 
tion in general (|Sion et al.lll988t iPauli et al.ll2006D . We 
conclude that unlike the old systems LP 400-22 and 
NLTT 11748 that show peculiar kinematics, the system 
GALEX J1717+6757 is kinematically young. The cal- 
culation of additional evolutionary tracks covering the 
observed range of parameters, particularly between 0.17 
and 0.19 Mq, would benefit studies of ELM white dwarfs 
and their likely progenitors. 
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